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Abstract 



We study the importance of thermal fluctuations during the electroweak 
phase transition. We evaluate in detail the equilibrium number density of large 
amplitude subcritical fluctuations and discuss the importance of phase mixing 
to the dynamics of the phase transition. Our results show that, for realistic 
Higgs masses, the phase transition can be completed by the percolation of the 
true vacuum, induced by the presence of subcritical fluctuations. 
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I. INTRODUCTION 



The dynamics of the electroweak phase transition has received much attention in recent 
years j^. The main reason for this interest is that, within the context of big-bang cosmology, 
the standard electroweak model can, in principle, satisfy the three conditions obtained by 
Sakharov for the generation of the observed baryon asymmetry in the Universe [Q: (a) 
baryon number violation; (b) CP violation, and; (c) nonequilibrium dynamics. To date, 
most mechanisms invoked to satisfy the third condition make use of a first order phase 
transition. [For mechanisms based on cosmic strings see Ref. |^.] This is also true for 
extensions of the standard model, which are currently favored by most authors, due to 
difficulties in generating sufficient baryon number within the minimal standard model [Q. 

Based on estimates of the quantum-corrected effective potential, nonequilibrium con- 
ditions are generated by the motion of critical bubbles of the broken (true) phase, which 
are nucleated within the symmetric phase. For the several baryogenesis models which rely 
on extensions of the standard model, the parameter space is large enough to justify the 
assumption of a fairly strong first order transition. However, analyses based on the 1-loop 
finite temperature effective potential for the minimal standard model and its improved 
versions, show that the phase transition is very weak, within the current lower bounds 
for the Higgs mass, ttih ^ 60 GeV or even of second order for larger Higgs masses, 
ttlh ~ 90 GeV P|. The weakness of the transition is further supported by nonperturbative 
lattice computations [|10|. 

Our interest in the present paper is to further investigate P,|TT| the possible consequences 
of having a weak first order phase transition at the electroweak scale. Although we will 
restrict our analysis to the standard electroweak model, our results can be adapted to any 
of its extensions. In fact, we will show that the strength of the transition can be used as 
a new constraint on the parameters of the model, always a welcome addition to the often 
large parameter space of extensions to the minimal standard model. 

The interest in exploiting the dynamics of weak first order transitions goes beyond its 
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potential relevance for baryogenesis. As recent results have shown, a sufficiently weak first 
order transition will exhibit a different dynamics from the usual homogeneous nucleation 



results [0; the extra free energy available in large amplitude fluctuations which are not 
included in the Gaussian computation for the nucleation rate will act to decrease the decay 
barrier, suppressing supercooling, and speeding up the completion of the phase transition. 
For even weaker transitions, critical bubble nucleation may be completely absent. Clearly, 
the information from the effective potential is not sufficient to determine the details of the 
transition; mean-field theory breaks down in the presence of large infrared corrections. 

In order to quantify the above statements, we will make use of the subcritical bubbles 
method P,pJ|. That is, we will model large amplitude thermal fluctuations by Gaussian- 
shaped bubbles of approximately correlation volume. Previous results based on a kinetic 
approach, have indicated that such fluctuations can destroy the ffist-order character of the 
transition for Higgs masses of order mn ~ 55 GeV Here we would like to complement 
this calculation by computing in detail the nucleation rate for such configurations, which 
was previously assumed for simplicity to be F = Am^{T)exp[—B/T], where A is a constant 
of order unity, m'^{T) is the curvature of the potential in the symmetric phase, and B 
is the free energy of the Gaussian configuration. Note that this is a nontrivial exercise, 
as these configurations are not solutions to the Euclidean equations of motion. Within 
reasonable approximations, we will be able to obtain the equilibrium number density of 
these configurations as a function of the tree-level Higgs mass, to show how the weakness of 
the transition is closely related to the breakdown of the dilute gas approximation. 

The remainder of this paper is organized as follows. In Sec. II, we give a brief review 
of subcritical bubbles in the context of the electroweak standard model. In Sec. HI, we 
compute the equilibrium number density of subcritical fluctuations and discuss the choice 
of parameters in the subcritical bubble configuration. There we also discuss the validity 
of the approximations taken and the range of applicability of our results. In Sec. IV, we 
discuss how the weakness of the electroweak phase transition is related to the breakdown of 
the dilute gas approximation; for large enough Higgs masses, subcritical bubbles percolate. 



completing the transition. In Sec. V, we have our main conclusions. Some technical details 
are left for an Appendix. 



II. SUBCRITICAL BUBBLES 



Following the work of Ref . |Tl|] , large amplitude fluctuations describing thermal fluctua- 
tions are parameterized as, 
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'Pscir) = <^A{T)exp y-^^j , (2-1) 

where (psc{r) describes (spherically symmetric) fluctuations in the scalar fleld, with amplitude 
(fA and radius given by R{T). The minimum value for R{T) should be compatible with the 
coarse-graining scale of the model. Later, in section III, we will come back to the parameters 
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O) and discuss our choice for (p^ and R. The fluctuations described by (|2.1|) approximate 



rather well the relevant fleld conflgurations, as recent work has shown [|13| 



A. Free Energy for Electroweak Subcritical Bubbles 

Let us estimate the free energy associated with the conflgurations given by We 
shall use, as a particular case, the electroweak model. 

In the electroweak model, for a Higgs self-coupling X <^ g'^ (with g denoting a generic 
gauge couphng), contributions due to scalar loops to the flnite temperature effective po- 
tential are small compared to the contribution due to gauge flelds and fermions. This is a 
common approximation employed in the literature which results in a 1-loop flnite tempera- 
ture effective potential given by [1^ 

1^(0, T) = D{T^ - T|)02 - ET03 + ^ ^2.2) 

where D and E are constants given in terms of the W and Z boson masses and of the top 
quark mass as: 
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and 



E 
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cr 



0.169 



(2.3) 
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(2.4) 



where a ^ 246 GeV is the vacuum expectation value of the Higgs field, mw = 80.6 GeV, 



mz = 91.2 GeV and we use rrit ~ 174GeV [|l^,|]. T2 in V{(f),T) (the spinodal instability 
temperature) is given by 
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where mj^ = (2A + 12_B)cr^ is the physical Higgs mass and B = g^^-^(6m|j/ + 3m| — 12m 



-0.00456. At in V{(j), T) is the efi^ective Higgs self-coupling (at 1-loop) given by 



At — A 
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(2.6) 



where the sums are performed over bosons and fermions, with degrees of freedom gi, and gj, 
respectively. In (^.61), Inc;, = 5.41 and Inc/ = 2.64. 

For temperatures T < Ti, where Ti is given by the solution of E'^T^ = |D (Tf — T2) Xti, 
V{(f), T) has minima at <y9 = and at 
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At the critical temperature T^, 
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ipf = describes the metastable 



we have V{(f) = 0,Tc) = V{(f) = v9+,Tc) and below Tc, ( 
phase (the false vacuum), while (/) = (p+{T) is the stable phase (the true vacuum). 

An important effect from higher loop corrections to V{(j),T) in (|2.2| ) is the reduction 
of the coefficient of the cubic term 0^ by 2/3 ||TB|. However, we are mostly interested in 
a toy model computation of the number density of subcritical fiuctuations, and will adopt 



the simplest 1-loop potential. For recent considerations on improving the 1-loop effective 
potential by resumming the most important infrared contributions and higher order graphs 
see, for instance, Refs. P jr7|JT8[] . 

The free energy for a given field configuration (psc{T) is, 
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F{T) = / d^x 



From (|2.1| ) and (|2.2|) , we get for F{T) the expression 



where 



F(T) = a{^A)R{T) + (3{^a)R\T) 
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D{r - Ti)^\{T) 
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III. NUMBER DENSITY OF SUBCRITICAL BUBBLES 



A. Partition Function for the Scalar Field in the Electroweak Model 

Let us define in the Electroweak model the partition function 



(3.1) 



where Xi denotes gauge and fermions fields (and ghost fields) and (j) is the SU(2) doublet 

(3.2) 
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where 0i {i = 1, . . . , 4) are real scalar fields. The tree level potential for the complex scalar 
field 0, given by 



yo(|0|) = V0V + A(0W, (3.3) 

for /i^ > 0, acquires a nonvanishing vacuum expectation value = a, which one assumes 
real and along, for example, the real component 03 of 0. Thus, in the broken phase, we 
define 03 = 03 + cr and 0i, 02 and 04 are the three Goldstone bosons. 

Let us denote by gi the coupling of the field with the Xi fields. If A ^ gf, i.e., the 
interactions among the field are weak compared with the — Xi interactions, then we may 
formally integrate out the x fields in to obtain 



D0e-^('^) , (3.4) 
where 

Vr(0) = -Iny D;^^e-/o'''^"/'^'"^E:uci.(*,x.) _ (3 5) 

For vector fields, the integration measure above includes the gauge fixing and ghost terms. 
We choose to work in the Landau gauge, which is the one usually used in the studies of the 
electroweak phase transition. Expanding W{(f)) in a derivative expansion, 

w{<j>)= f^dr [d'x\vom+vm)+zm{d,<py{dy) + ...] , (3.6) 





where Vb(|0|) is the tree level potential ( p.3|) and Va(|0|) is the contribution of the Xi loops, 
coming from the integration over the Xi fields in with the scalar field in the external 

legs. Since V{(f), T), Eq. ( p.2|) , is obtained by neglecting scalar-boson contributions, which is 
analogous of just making a functional integration on vector and fermion fields at the 1-loop 
approximation, Vq + Vg above, at the 1-loop approximation for the Xi fields, can be written 



as in Eq. ( p.2|) . 2'(|0|) is the wave-function renormalization factor. 2'(|0|) has already been 
evaluated by many authors, in the determination of an effective action for the Higgs field 
with the objective to study corrections to critical bubbles nucleation rates [see, for instance 
||T9|-pTl], where it is utilized a process of integration of fields, or field degrees of freedom, 
analogous to the one done above to determine W{(j)). Although 2'(|0|) receives nonvanishing 



1-loop contributions, it has also been shown in [|I9| that these contributions are expected to 



yield only small corrections to the effective action and, therefore, these corrections could be 
treated as perturbations. Later, in section III.D, we show that the same approximation can 
be taken here, at least for the range of temperatures and higgs masses we are interested in. 
Therefore, as a first approximation, we will neglect all wave function corrections and take 
Z(|0|) ~ 1, for simplicity. Thus, at 1-loop order, we write the "effective" action W{(j)) as 
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W,{(P)- I dr I d'x \d^^\'' + V{\(PlT) . (3.7) 







B. Number Density for Subcritical Fluctuations 

Well-known results show that, in a dilute gas approximation, the average number of 
extended objects (for example, topological defects) described by some field configuration cpc 
can be given by |2^-^ 

where Z{ipc) [Z{fv)] is the partition function of the system computed by expanding the 
scalar field around the field (vacuum) configuration ip^ {<Pv)- For stable configurations the 
ratio in (|3.8|) is real. However, for unstable configurations (like, for example the sphaleron. 



the critical bubble or bounce configuration and also for ipsc, given by ( pTTP ), the ratio in ( p.8|) 
is complex due to the existence of negative eigenvalues, associated with the instability of the 
configuration. This is the case for subcritical fluctuations. In this case, we will adopt the 
procedure of Arnold and McLerran in [^] for the case of sphaleron configurations, where 
they associated the average number of sphalerons to the total rate of transitions multiplied 
by the time of a single transition, giving, in the dilute gas approximation, 

A'. ~ ~ h.^^ , (3.9) 

where V is the transition rate given by F ~ ^Im^|^ and uj^ is the negative eigenvalue. 
Note that in this case the final result for ( |3.9| ) will not depend on the negative eigenvalue, 
as shown in Taking ( |3.9|) as also valid for the case of subcritical fluctuations, we can 



associate T to the nucleation rate of subcritical fluctuations and the possible negative eigen- 
values, which must appear in (p.8|), are the ones associated with the collapse mode of the 
fluctuation. We will, therefore, evaluate Nc for subcritical fluctuations not taking into ac- 
count the possible imaginary eigenvalues associated with the instability of the conflguration, 
that is, we will adopt a procedure similar to the one expressed by ( p. 91 ). 

Computing the partition functions will give us the equilibrium values for relevant physical 
quantities. In particular, the equilibrium number density density of subcritical fluctuations, 
we will thus write as 

where rj = ri{x,T) and ^ = ^(x, r) are small perturbations around the conflgurations ipsc 
and (which we take as the false vacuum conflguration (ff = 0), respectively. In the 
following section we will discuss the limits of applicability of Eq. ( p.lO| ) within the standard 
electroweak model and in other situations applicable to phase transitions in general. 



C. 1-loop Evaluation of A's, 



In the 1-loop approximation for W{(f)) in (|3.10|) , given by (|3.7|) , Wi{(f>) can be expanded 
about a fleld conflguration ip^ as 



+ ^ / d\d^x'W';{^c] X, x')7]{x)7]{x') + Oir]'') , (3.11) 



2! 

where ri{x) = (f){x) — ^{x) and / d'^x = Jq dr J d^x. Taking ri{x) as small perturbations 
around the fleld conflgurations ip{x), the terms of order 0{ri^) and higher can be treated as 
small perturbations. Note that for the kind of conflgurations we are dealing with, W[{ip) = 
'^^0 I 't>='P "io^s not vanish in general {Lpsc{r) since is not a stationary solution of Wi{(l)). 
Using ( 3.11 ) in ( ^.l(j| ) we can perform the gaussian functional integrals and since we 



have already integrated out all other flelds interacting with the scalar field (p, the scalar field 
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propagators will include loop (quantum) corrections from the other fields coupled to 0. We 
must, therefore, take some care when performing the gaussian functional integral in in 
order to avoid possible double-counting. From ( |3.11| ) and ( p.lO| ), with ipy = ipf = and 
Wl{(py) = 0, we get 



det Wi'iif, 



(3.12) 



where AWi = Wi{ipsc) — Wi{ipv) and W{' denotes the correct dressed (including the 1-loop 
quantum corrections from the Xi fields in (|3.1| )) inverse propagator for the (p field: 



(3.13) 



with m|(v9) = 2D(T'^ —T2) +3Xt^p'^ ■ Note that mp{({>) is the classical Higgs mass corrected by 
the self-energy (T 7^ 0) corrections coming from the Xi fields (fermions and gauge bosons). 
mp{i~p) does not include the term proportional to E of ( p.2|) , which would give origin to 
a linear term in ( p.l3|) . This term is absent since it would not appear in the self-energy 
corrections to the scalar field and also that the presence of such a term in ( p.l3| ) is well 
known to lead to a wrong counting of loop corrections to the scalar field effective potential 
(see for instance, ref. [jT^)- 
In ( ^.121 ) we also have that 



d''xd''x'W[{ip,,-x)W[{ip,,-x'){x\ Wl{ip,,-x,x') \x') 



(3.14) 



where 



{x\ Wl{^^^]x,x'] 



-1 1 +~ 

n=-c 



^3j^ ^iuJn(T-T')+ik.{x-x') 



(3.15) 



are the Matsubara frequencies. In appendix A we solve 



where w„ = ^ (n = 0, ±1, ±2, . 
( CT ) exphcitly. 

In order to compute the determinant ratio in ( |3.12| ), we must first isolate the possible 
zero modes in det Wi{ipsc) = det[— □ + m|(y9sc)]- Note that, as exposed above, here we will 
not take into account possible imaginary eigenvalues. Noting that iyi(0), Eq. ( p.7|) , and the 
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corresponding field equation, = 0, are translational invariant. Therefore, even tliougli if 
V^sc is not a solution of the field equation, there must be three zero eigenvalues associate to 
(fsc, related to the three translational modes. We also have three more zero modes that are 
related to the rotational symmetry SU{2), associated with the three Goldstone bosons in the 
broken phase. We handle these zero modes by the standard way, by introducing collective 
coordinates, such that the determinantal ratio in ( p.l2| ) can be written as: 



det Wi'i^,) 



trans. 'rot. 



det'Wj'jifs 
det Wi'{ip, 



(3.16) 



where fitrans. is the usual factor coming from the translational modes, given by 



trans. 



271 



V 



(3.17) 



and firot; due to the rotational modes, has been explicitly obtained in for the case 
of fluctuations around the critical bubble configuration, and in our case we can write the 
analogous expression for the (psc configuration as 



O - — 



(3.18) 



In ( p.l6|) the prime in the determinant is to indicate that the six zero eigenvalues have been 
excluded. In (|3.17|) V is the space volume. 

From ( p.7|) , we can write the determinantal ratio in ( p.l2| ), in momentum space, as 



exp< --In 



(3.19) 



where the productories in i and j are, formally, over the eigenvalues Ef{ipsc) and i?|((y9^), 
respectively. In ( |3.19| ), we have used the identity In det O = trlnO. Taking into account the 
six zeroes modes of ( p.l9| ), associated with the three Goldstone bosons of the Higgs doublet, 
in the broken phase, and the three translational modes, we obtain 



det lyf (<^s, 
det W^{^^ 



^^trans.^^rot. GXp {-Ij^^^^ 



(3.20) 
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where the factors f^trans. and Qrot. are given by ( p.l7|) and ( p.l8| ), respectively. The prime 
in Yli in ( |3.20| ) is a reminder that those zero modes, associated with (p.l7| ) and ( |3.18[ ) have 
been excluded from the product of eigenvalues. 
Using the identity: 



+00 



n 1+ 



n=l 



(27m) 



sinh(a/2) 
a/2 



(3.21) 



we get for (|3.2CI|) the expression 
"detWn^" 



' j 1 

~ ^trans.^rot. CXp i — 



-i3Ei{if sc) 



2E 



|E,(v..)+ln(l-e-^^^-(-'') 



(3.22) 



For the field configuration (p^c, given by ( |2.1|) , and the vacuum configuration (py (the false 
vacuum, = 0), from (p.7|), we can write that 



E 

(continuum) - 



(^Ei^) + In (1 - e- 



(fx 



(2^ 



+ In 1 - e 



(3.23) 



where m|((/?) = 2D{T'^ - T|) + 3At</?^. Note that the terms like Jd^kyk^ + m? in (|]23D, 
that are ultraviolet divergent, can be subtracted by introducing the usual counterterms of 
renormalization for the scalar field loops, rendering the exponent in ( p.22|) finite. 
Using (|3.23|) in ( |3.22| ), Eq. (|3.12|) can therefore be written as 



Z(y(psc) 1— loop approx. 
Z{(fv) 



^trans.^rot.^ CXp 



AFcff(T) 
T 



(3.24) 



where AFefr(T) denotes an effective free energy for subcritical bubbles, given by 



AFeff(T) = F(T) + j d^x [V^ivscT) - \/^(<^.,T)] 



- T / d^xd^x'Wii^sc, x)W[{^sc; x){x\ W['{^sc] X, x') \x') 



(3.25) 



where F{T) is given by ( |2.10| ). The scalar field quantum contribution in ( p. 251 ), in the high 
T limit, is given by 
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'-Xt^\{T)T^R\T) 



32 



(3.26) 



The last term in ( |3.25D , from ( |3.14D and ( |3.15| ), can be written as 



/ 



-vVsc(x) + \/'(</^sc(x),r) 



-V' ^,,{x') + V'{^,,{x'),T) 



J(|f-f|), (3.27) 



where I{\x — x'\), in the high temperature hmit, is given by (see the Appendix) 



l[\x — x\) ~ 



TT 



ma x—x 



(3.28) 



(27r)2|f-f' |T 

Using the above expression in ( |3.27| ), we can see that the largest contribution will come 
for values of \x — x'\ close to the inverse of rrijj, i.e., for values close to the correlation length 
^(T) (~ l/m/j). We therefore may restrict \x — x'\ to the size of the subcritical fluctuations 
and we here consider R{T) ~ C,{T)), thus obtaining 



I{\x-x'\) 



IT 



{2'KfTR{T)e ■ 

Substituting (y5sc(x) in ( |3.27D and using the above expression for I{\x — x'|), we get 



(3.29) 



^ ' Te 



D{T' - T|) - 3^ET^a{T) + ^Xt^UT) 



. (3.30) 

Z{<psc) 



8 ^ ' 18 

Using the expression for F(T), Eq. (|2.1(]| ), and Eqs. (|3.26|) and (|3.30|) , we get for ^ 
the high temperature limit and at 1-loop order, the expression 



m 



VT-'A{T)exp[~B{T)] 



(3.31) 



where A{T) and B{T) are given by 
A{T) 



'-'~e*"fi»(T)^«(T) 
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^ + y^DiT' - Ti)R\T) - 1^ETMT)R\T) + 
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+ ^vliT)R\T) 



(3.32) 



and 
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B{T) 



l-loop, highT AFcfj(T) 

T 



(3.33) 



with AFcfr(T) given by ( ^1251) together with Eqs. ( ^ ), j ^M ) and (lOOl) . 

From ( p.8|) , the equihbrium number density of subcritical fluctuations, in the dilute gas 
approximation, is given by 

1 Z{(psc) 



risciT) 



V Z{^,) 



(3.34) 



D. A Discussion of the Validity of the Approximations and the Choice for the 

Parameters and R 

We expect that the expression for ngc, obtained in the dilute gas approximation, be valid 
as long as the following holds 

[nsc(T)]^i?(T) <1, (3.35) 

We can also express the condition in terms of the volume occupied by the fluctuations 
described by ipsc, at some flxed temperature T, Vsc{T), 

Att 

Kc(T) = —R\T)n,,{T)V , (3.36) 



such that (|3.35|) can be reexpressed as 



< 1 , (3.37) 



V 

that is, the dilute gas approximation is a good approximation as long as the volume occupied 
by the subcritical fluctuations be small. 

Another approximation that we considered, the neglected of wave function corrections 
to the effective action, can be thought to be somewhat problematic. However, we can 
estimate their relative contribution to our results and inquiry whether our approximation of 
neglecting these contributions is good enough. For such estimation we may take the leading 
order correction to AFcff(T), due to the wave function correction coming from the functional 
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integration over the gauge boson fields. As an order of estimative, we may take tlie result 



obtained in |T9[, for the SU(2)-Higgs model, which we write as, in neglecting contributions 



due to plasma masses. 



Z{^,T) = ^^—, (3.38) 



where g = 2my^ /a ~ 0.328. An estimate of the relative contribution of the above factor can 
be given by the ratio (as in |1T9| ) 



sc ) llg T 

= 2 ~ — /= 1 (3.39j 

/ (Px (V<^sc) 4v 27r 



where we have used ( p.l|) in the rhs of the above equation. It is interesting to note that the 
above expression do not depend on the subcritical bubble configuration radius R. 

The range of applicability of the dilute gas approximation and the validity of the ap- 
proximation of neglecting wave function contributions is determined once the parameters 
if A and -R, of our ansatz, Eq. ( |2.1| ), are set. 

In the following, we take the amplitude as been given by the temperature dependent 
broken minimum, ip^{T), Eq. ( |2.7| ). For R we assume it as been given by the correlation 
length ^(T), given by the inverse of the temperature dependent mass in V{(f),T). We next 
justify these assumptions. We expect that fluctuations in the scalar field would appear in the 
system with arbitrary amplitudes and sizes. However, physically, we can set some limits on 
these parameters. For example, fluctuations with too small amplitudes are already summed 
over in the computation of the (coarse-grained) effective potential. The relevant fluctuations 



for the dynamics of the phase transition can roughly be identified [p!2tp!3| as been those with 
amplitudes ipA ~ fmax, where v^max is the value of the higgs field at the maximum of the 
effective potential. Among these fiuctuations in the false vacuum, those with (pA close to 
the local minimum cp^ are certainly the most expected. About the fluctuations radius, R, 
their minimum radius must be compatible with the coarse-graining scale of the model and 
the coarse-grained effective potential. In general, we can expect that R > C,(T). Here, 
we implicitly assume that the most probable fluctuations in the system are those with 
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radius close to the correlation length ^(T). A better approximation would be taking R as 
the average radius of the fluctuations. Unfortunately, a reliable method to determine this 
quantity is still lacking. Recently, the authors of ref. ||2^ propose a method to determine 
the average radius R of fluctuations, considering i? as a truly dynamical variable, however, 
the obtained result for R is much smaller than ^(T). In another recent method proposed by 
the authors in [|l2|,|l3l, a statistical mechanical method is used for studying the importance 



of phase mixing in the phase transition, where the volume fraction occupied by fluctuations 
is determined by summing over all possible fluctuations of different sizes and amplitudes. In 
the next section, by computing this volume fraction within our approximations, we will be 
able to compare our results with the ones obtained in the statistical mechanics method. 
With (fA = f+{T) and R = ^(T), we can immediately evaluate (|3.39| ) and have an 



estimative of the leading wave function contribution. In Figure 1 we have given 5 as a 
function of the higgs mass, with 6 computed at the critical temperature Tc. Our result for 
the fraction given by ( |3.39| ) is much like the same one as given in |T^, where S > 1 for 
^ SOGeV. However, for such larger higgs masses it is well known that the perturbative 
expansion breaks down and the effective potential becomes unreliable [|7|,PJ2^. For the 
interval of higgs masses we will be interested in, the wave function correction is small and, if 
higher order corrections are taken into account and plasma masses are included in Z{(j),T), 



we expect the wave function contributions be even smaller, as shown in |19 



IV. PERCOLATION OF SUBCRITICAL FLUCTUATIONS 

Studies of the electroweak phase transition, through the analyses of the effective potential 
26| and direct numerical simulations p8[, indicate that the nature of the electroweak phase 



transition is most possible of first order, but it has a excitation barrier to small, qualifying 
the phase transition as very weakly first order. In these situations it has been argued that 
the dynamics of the phase transition could be quite different from the usual process of 
critical bubble nucleation. In fact, due the weak nature of the phase transition, (subcritical) 
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fluctuations over the barrier can become the main mechanism by which the phase transition 
can complete. Here, we analyze the importance of a possible large amount of phase mixing 
at temperatures higher than the nucleation temperature of critical bubbles. In this sense, 
we argue that if the volume of true vacuum phase fluctuations becomes higher than a certain 
value, percolation of the fluctuations can occur. We thus have a possible scenario where the 
phase transition can be completed just by those over the barrier subcritical fluctuations. 

A well known result of statistical mechanics of dynamics of cluster systems |^9| show 
that there is a critical probability value pc, with the cluster probability defined by 

Cluster Volume , , ^ , 

P=^, ^71 ' ^^-^ 

bystem Volume 

where, in three space dimensions, this critical percolation probability is roughly ~ 0.3 
29| and beyond which clusters are favored to coalesce and grow, filling the whole volume of 



the system. 

From our definition of the volume occupied by the field configurations (psc, at a fixed 



temperature T, Eq. ( |3.36| ), it is quite natural to define, therefore, a percolation temperature 
Tp as given by 

Kc(Tp) 



V 



0.3. (4.2) 



From Eq. ( |3.36|) , we obtain 



y/?3(TpKc(Tp)~0.3. (4.3) 

Condition (|4.3| ) gives us the temperature (as a function of the Higgs mass) below which 
subcritical fluctuations become favorable to coalesce and grow, forming large regions of the 
true vacuum phase. If Tp is larger than the nucleation temperature T/v, that marks the onset 
of critical bubbles nucleation, then we have an effective mechanism by which the phase tran- 
sition can complete, as explained above, only due the dynamics of subcritical fluctuations 
of the stable phase ip+, inside the false vacuum phase. Note that at Tp, condition ([4.2[), 
due to (|3.37|) , will also signal the breakdown of the dilute gas approximation. Thus, for 

16 



temperatures too close to Tp, the obtained expression for Usc, Eq. ( p.34|) , would just give a 



qualitative indication of the importance of subcritical fluctuations during the phase transi- 



tion. We must, however, comment on the alternative method of the authors of refs. |T2| , p!3 
where a statistical mechanical approach is used for studying the dynamics of subcritical bub- 
bles, modeled there as in Eq. ( |2.1| ). By constructing a kinetic equation for the (dynamical) 
number density of fluctuations, they could obtain an analytical expression for Vsc/V, within 
appropriate approximations. Their kinetic approach also break down for Vsc/V close to the 
critical percolation probability, where fluctuations become more dense and the treatment 
of the coalescence of fluctuations becomes important. Besides the different approach, the 
method applied to the electroweak phase transition (see Gleiser's talk in |jl3[) show that 
Vsc/V changes sharply for ttih ^ GOGeV, just as the result obtained here for Vsc/V shows 



(Figure 2). Also, recent numerical simulations of very weak first order phase transitions ||28 
give support to our results. All these dynamical studies indicate the validity of the various 
approximations here adopted in order to make possible to arrive to an analytical expression 
for Use- 

In our analyses it is enough to study the ratio Vsc/V at the critical temperature T^, for 
which we get the plot shown in Figure 2, in terms of the Higgs mass mn, where we have 
used the approximation for At, 

From Figure 2 we see that for a Higgs mass of ~ 60GeV, our results indicate that phase 
mixing is so large that the volume of subcritical fluctuations in the broken phase start 
becoming relevant for the dynamics of the phase transition. However, due to the condition 
of validity of the dilute gas approximation used in our evaluation, Eq. ( |3.37| ), our results can 
only be interpreted qualitatively for itih ^ GOGeV. Besides these limitations, the results 
obtained here are in accordance to recent results based on dynamical studies of very weak 
phase transitions, as discussed above. 

The main result of this paper is the computation of the equilibrium density for the field 

17 



configuration (psc, with the full evaluation of the preexponential term A{T) given by (p.32|) 
and the effective free energy AF{T), Eq. ( |3.25| ). Previous results have given just in 
terms of a Boltzman distribution form, with free energy given by ( p.lO| ) and prefactor given 
by just T^. From ( |3.31| ) and ( |3.34| ) we see that the preexponential factors differ by the factor 
A{T), which we show in Figure 3, as a function of the Higgs mass, which shows that the 
prefactor can assume large values for small values of the higgs mass. In Figure 4, we also 
compare the exponential factor AFcs{T) with F{T), both compute at T = Tc. 



V. CONCLUSIONS 

We have explicitly evaluated the equilibrium density of subcritical fluctuations and shown 
that the same can differ substantially from the usual expression. We have also analyzed 
the relative importance of subcritical fluctuations during the electroweak phase transition, 
showing that the mechanism for completion of the phase transition can be quite different 
from the usual critical bubble nucleation in a first order phase transition. For values of the 
Higgs mass at and higher than the experimental lower bound of ~ 60GeV, our results show 
that substantial phase mixing can be present prior to critical bubble nucleation, changing 
the dynamics of the phase transition. Our results for the number density and volume 
occupied by those subcritical fluctuations can be interpreted as an average value at some 
fixed temperature, since we expect that in some large volume V ^ R^(T) (-R(T) ~ C,(T)) we 
have both processes of nucleation of subcritical bubbles and their shrinking, keeping their 
average density, at that temperature, constant in that large volume. 

We must also remember that the results obtained here are in accordance with other stud- 
ies involving the relevance of thermal fluctuations during the electroweak phase transition 
1^ and recent studies of the dynamics of very weak first order phase transitions, both by nu- 
merical simulations as by a statistical mechanics approach for the dynamics of fluctuations 
modeled by ipsc, Eq. (|2.1|). 
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APPENDIX A: EVALUATION OF EQ. (3.28) 



In this appendix we evaluate (|3.14|) and get the expression for I{\x — x'\), Eq. ( p. 281) . 
The term expressed in ( p.l4|) , 



W^scM W^sc= d^xd^x'W[{y^,,;x)W[{^,,;x'){x\ W^{v,,-x,x') \x') , (Al) 



where W[ = is obtained from Eq. ( |3.7] ). Since W[{(f)) computed at = ^Psc{x) is 
independent of the Euchdean time, the two Euchdean time integrations in apphed 
directly to (x| Wi{(psc'-, x^x') \x') which we called I{\x — x'\): 



[P [1^ , \ - I, -1-1 
Jo Jo ^ J 



(A2) 



Using ( |3.15|) , performing the sum in the Matsubara frequencies in ( |3.15|) and using the result 
in we can perform the two Euclidean time integrations in (|A2|) in order to get the final 
result 

d^k Qik.{x-x') 



I{\x-x'\) 



2 P 



m 



3 

/3 



sinh f3\ k"^ + m 



'/3 



(A3) 



Performing the angular integrations and the k integration of ([A3|) , we get, in the high 
temperature limit, Pm^ <^ 1, the result given in Eq. ( p.29|) . 
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Figure 1 
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Figure 1: An estimate of tlie relative magnitude of the leading wave function contribution, 
Eq. ( |3.39| ), computed at T = Tc. 
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Figure 2 
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Figure 2: The volume fraction ^ of true vacuum fluctuations (computed at the critical 
temperature Tc) as a function of the Higgs mass mn (GeV). 
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Figure 3 
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Figure 3: The Prefactor coefficient A{T), computed at the critical temperature Tc, as a 
function of the Higgs mass. 
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Figure 4 
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Figure 4: The exponential factor B(T) = ^^''^^'^^ (upper curve) compared with F[T)/T 
(lower curve), both evaluated at T^,, as a function of the Higgs mass. 
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